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In this paper, the time series analysis method CLEANest is employed to search for characteristic periodicities in the radio light 
curves of the blazar NRAO 530 at 4.8, 8.0 and 14.5 GHz over a time baseline of three decades. Two prominent periodicities on 
time scales of ~6.3 and ~9.5 yr are identified at all three frequencies, in agreement with previous results derived from different 
numerical techniques, confirming the multiplicity of the periodicities in NRAO 530. In addition to these two significant periods, 
there is also evidence of shorter-timescale periodicities of ~5.0 yr, ~4.2 yr, ~3.4 yr and ~2.8 yr showing lower amplitude in the 
periodograms. The physical mechanisms responsible for the radio quasi-periodic oscillations and the multiplicity of the periods are 
discussed. 
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1 Introduction 



The energy production and release of active galactic nuclei 
(AGN), which represent a population of the most powerful 
sources in the Universe, is one of the hottest problems in ex- 
tragalactic astrophysical studies. The vast energy supply of 
the AGN is thought to be related to the accretion of the cen- 
tral supermassive black holes (SMBHs). Imaging of the en- 
ergy engine of the AGN, i.e., the BH-accretion disk system, 
which is extremely compact, is still limited by the angular 
resolution of the current observing facilities. Due to the diffi- 
culties of the imaging observations, non-imaging techniques 
are extensively used for investigating the dynamic processes 
taken place in the galactic nucleus. For instance, the flux den- 
sity variability of AGN, which is rather easy to acquire from 
monitoring observations with a single telescope, provides vi- 
tal information to constrain the energetic processes related 
to the black hole system. Moreover the time series analysis 
of the variability data is done in time or frequency domain, 
and thus it is free of angular resolution limitation. Blazars 
consist of a sub-class of AGN and are characterized by rapid 
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and violent variability across the whole electromagnetic spec- 
trum. Therefore variability studies of blazars open a window 
to understanding the AGN physics. Searching for periodicity 
from the variability of blazars is of particular interest, be- 
cause the inferred time scale may serve to constrain the size 
and location of emitting source and the dynamic parameters 
of the physical models of AGN (e.g., [1]). Variability analy- 
sis of blazars from radio, optical, X-ray bands on diverse time 
scales has been reported in literature (reviewed by [2-6]). 

Many numerical methods that were originally developed 
in statistics, signal processing and communication sciences 
have been applied to study the periodic oscillations in astro- 
nomical light curves. A basic tool of periodicity analysis is 
the Discrete Fourier Transform (DFT). However, the ordinary 
Fourier transform may often induce fake peaks in the power 
spectrum when unevenly sampled data dealt with [7-10]. Due 
to various observing problems, light curves of astronomical 
objects can be regarded as irregularly sampled time series. 
When the astronomical time series have large gaps, that of- 
ten results in sharp spikes in the periodograms. The noise 
of the sparsely sampled data may also result in fake peaks 
in periodograms. In addition, the observed data may contain 
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multiple modes of periodic signals, which interact with each 
other and make the periodogram even more complex. An al- 
gorithm named CLEANest [11] has been designed on the ba- 
sis of date-compensated discrete Fourier transform (DCDFT: 
[12]), and it is proven to be a robust and efficient technique 
to identify the real periodicities from noisy data containing 
multiple periodic signals (e.g., [11, 13]). The DCDFT and 
CLEANest method have been used for period detection of 
blazars [14-19] from the optical and radio light curves. The 
derived periods range from 2.2 to 20.8 years [15]. 




Figure 1: The radio light curves of NRAO 530 at 4.8 (square), 8.0 (cir- 
cle) and 14.5 GHz (triangle). In order to avoid the overlapping of three light 
curves, the 4.8-GHz flux densities have been subtracted by 2 Jy, and the 14.5- 
GHz flux densities have been added with 2 Jy. The strongest outburst occurs 
during 1996-1997. 



radio telescope of the University of Michigan Radio Obser- 
vatory (UMRAO). The observations were made at three fre- 
quencies, i.e. 14.5, 8.0 and 4.8 GHz. The light curves are 
exhibited in Figure Q] The starting epoch of the monitor- 
ing observation is April 1, 1978 (4.8 GHz), April 19, 1967 
(8.0 GHz) and October 10, 1974 (14.5 GHz). And the latest 
epoch is in November 2009 for all three frequencies. The to- 
tal time span is 31.6 (4.8 GHz), 42.7 (8.0 GHz) and 35.0 yr 
(14.5 GHz), which is much longer than the radio data used 
in [15] by 50%, 33% and 46%, respectively. Accordingly the 
data points are more than [15] by 22% (4.8 GHz), 10% (8.0 
GHz), and 16% (14.5 GHz). The observations before 2001 
were regularly performed; an average gap between the sam- 
pling is about 2 weeks. After 2001, there were two big gaps 
appearing in 2001-2003 and 2005-2008. Considering that 
the typical periods of the radio variability of blazars are of 
the order of a few years and the strongest outburst lasts a few 
months, the biweekly sampling is good enough for perform- 
ing an accurate periodicity analysis. 

This paper will make use of the CLEANest method 
to search for periodicities from the radio light curves 
of NRAO 530. Benefitting from longer time series and 
more data point, the present work aims to detect possible 
weaker period(s) and to confirm the multiplicity of the pe- 
riodicities. The structure of the paper is organized as follows. 
In Section 2 the CLEANest algorithm is briefly described. 
Section 3 presents the results from the CLEANest method, 
and interpretations of the observations are given in Section 4. 



NRAO 530 (PKS 1730-130) is a well-known blazar at a 
redshift of 0.902 [20]. The light curves of NRAO 530 show 
prominent outbursts across the whole electromagnetic spec- 
trum, and VLBI observations reveal apparent superluminal 
motions of the jet knots, indicative of relativistic beaming 
of the jet flow ([21] and references therein). The radio light 
curves display a number of outbursts, and the peaks show 
an apparent 6-yr cycle since 1980 (Fig. Q]). Previous vari- 
ability analysis of the radio data led to detection of two pe- 
riodicities of 10-yr and 6-yr (DCDFT: [15]; phase dispersion 
minimization: [22]). Except for these two most prominent 
periods, there were also hints of short-timescale periods [22], 
however the amplitudes of which were too low to be distin- 
guished from the noise. These shorter periods need further 
investigation to verify their reliability. Another interesting 
finding is the observations of multiple periods in NRAO 530. 
Although multiple periods have already been observed in a 
number of blazars [15, 19], it is still a puzzle whether the 
multi-period is a common feature of the blazar variability, 
and the physics of the multiplicity is not well understood yet. 
The blazar NRAO 530 deserves a comprehensive variabil- 
ity study owing to the above reasons and serves as a tem- 
plate to better understand blazar physics. In this paper, 
we shall study the periodic variability of NRAO 530 by using 
the monitoring data observed with the 26-meter paraboloid 



2 The CLEANest technique 

As mentioned in Section 1, classical Fourier transform of as- 
tronomical data which are irregularly sampled may often in- 
duce spurious peaks in the power spectrum due to large sam- 
pling gaps, confusing the identification of real periods. These 
fake peaks appear at frequencies oj s + ka> g , where cj s is the 
signal frequency, u> g is the frequency of the gaps, and k is an 
integer [11]. The noise of the observing data may also gen- 
erate fake peaks. Therefore discrimination of real periodic 
signals from noise is the primary task of period detection. 
Light curves of a substantial blazars contain multiple peri- 
ods, which are interacting and blending in a complex matter. 
In addition, some weak short-lived periodicities might be im- 
mersed in other strong persistent periodicities. Therefore a 
further consideration is to effectively extract individual peri- 
odic components. Although there are many existing numer- 
ical techniques that can be used for periodicity analysis of 
time series data, the CLEANest is found to be one of the ex- 
cellent techniques which is powerful not only at identifying 
noise-induced signals, but also in separating multiple compo- 
nents. 
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2.1 The principle of the CLEANest 




-10h "i i ; ; 

1975 1980 1985 1990 1995 2000 2005 2010 




Time (year) 



Figure 2: The simulated time series mixed with three periodic compo- 
nents of 10, 6, and 4 yr. The top panel shows the multi-periodic signal with- 
out noise and gaps, and the bottom panel shows the signal with noises and 
irregular gaps added. 

For periodicity analysis of unevenly spaced data, a widely 
used analysis technique is the modified periodogram [9, 23], 
which is based on a least squares regression onto the two 
trial functions (sin(wf), cos(wf))- The DCDFT is an improved 
Fourier-based periodogram, and is superior for dealing with 
the irregular sampling problems compared with ordinary Dis- 
crete Fourier Transform. It carries out least square regression 
fitting of the observing data with a linear combination of three 
trial functions, sin(wf), cos(wf) and a constant. The DCDFT 
improves the estimate of both the signal frequency and am- 
plitude. 

The CLEANest algorithm is exploited with an attempt to 
eliminate the fake peaks in periodograms resulting from ir- 
regularity of the time series. It is in principle a combination 
of the DCDFT algorithm [12, 1 1] and the CLEAN algorithm 
[24, 25]. An accurate identification of the signal frequency 
and the following-up subtraction of the Fourier component(s) 
from the spectrum is achieved by utilizing the CLEAN tech- 
nique [24]. In the CLEANest procedure, the strongest sig- 
nal peak together with the corresponding sidelobes are first 
subtracted from the original DCDFT spectrum, and a resid- 
ual spectrum is produced. Then the secondary component 
is searched for in the residual spectrum. If the remaining 
strongest peak in the residual spectrum is statistically signif- 
icant, the original data are analyzed again to find the pair of 
frequencies which give best fits to the data. These two Fourier 
components (primary and secondary) and corresponding fake 
components (associated with their sidelobes) are subtracted 
again, producing a secondary-level residual spectrum. This 
DCDFT-CLEAN cycle repeats until the residual spectrum is 
considerably flat and does not show any significant peaks. 
Since the frequency components involved would be refined at 



each run of the CLEAN procedure, under most circumstances 
higher-order CLEAN spectra would give progressively bet- 
ter frequency estimates. The power spectrum obtained in 
this way is called the CLEANest periodogram. The discrete 
CLEANed Fourier components are used to model the origi- 
nal dataset and create a reconstructed light curve which can 
be compared with the observed one. 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Frequency (yr 1 ) 

Figure 3 : The CLEANest periodograms of the multi-periodic simulated 
signal (shown in Fig. ff] bottom panel). From top to bottom, it displays 
the DCDFT, the first-iteration CLEAN, the second-iteration CLEAN and 
the final-iteration (i.e., the CLEANest) spectra, in sequence. Obviously, the 
CLEANest(3) successfully extracts all trial periodic components. 

2.2 Testing the CLEANest algorithm with a simulated 
multi-periodic time series 

In order to test the ability of the CLEANest in extracting the 
underlying periodic signals from the unevenly sampled time 
series, we applied the CLEANest method to the simulated 
multi-periodic signals of y — sin(j7Tt)+2 sm(^nt)+3 sin(|flt)> 
a combination of three sinusoidal functions with periods of 
4, 10, 6 yr. These three Fourier components have magnitude- 
less amplitudes of 1, 2, 3 respectively. The simulated light 
curve is shown in the top panel of Figure [2] The time span 
of the simulated time series is from 1974 to 2009, similar to 
that of NRAO 530. Standard normal distribution noise and 
random large time gaps were also added to the simulated sig- 
nals, as shown in the bottom panel of Figure [2] 

The DCDFT scanned the whole light curve and modeled 
the time series with two sinusoidal functions and a constant. 
The derived periodogram is shown in Figure [3] The high- 
est peak is located around a frequency of 0.167 yr 1 with 
an amplitude of ~2.5, in rough agreement with the strongest 
(amplitude = 3) trial period of 6.0 yr. To the left and right 
of the primary peak, there are two lower-amplitude (about 
1.5) peaks that are symmetrically distributed with respect to 
the primary peak. These two peaks are actually the sidelobes 
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of the strongest peak. The secondary peak in the DCDFT is 
at ~0.35 yr" 1 with an amplitude of ~ 1 .65. We run an itera- 
tion of CLEAN to subtract the primary peak from the original 
spectrum, and obtain the residual spectrum which is shown 
in Figure [3}b. The most significant component is at a fre- 
quency of 0.0957 yr" 1 and has an amplitude of 1.8. This 
Fourier component corresponds to the secondary strongest 
trial period of 10 yr. The amplitude of the peak at 0.35 yr" 1 
in Figure [3]-a is significantly reduced, suggesting that it is 
a high-order harmonic component of the strongest periodic 
component generated by the Fourier transform. Therefore 
it is not identified as a real periodic component. Moreover 
the sidelobes of the 0.167-yr" 1 primary component are also 
suppressed down to an amplitude level of 1.0. Repeating 
the CLEAN process, we successively obtain periodograms 
of CLEANest(2) and CLEANest(3). The residual spectrum 
of the CLEANest(3) becomes uniformly flat, having no sig- 
nificant spikes. CLEAN procedure stops at this step, and all 
three trial periods are distinctively identified. Note that the 
small shift of the detected periods away from the assumed 
ones reflects the intrinsic shortcoming of Fourier Transform 
when irregular noise is dealt with. We can see, after each 
run of the CLEAN procedure, the fake high-order harmonic 
component and the sidelobes of the real signal are effectively 
extracted from the original spectrum, and they would have no 
much effect in the residual spectrum. 



the higher-order harmonics of the 9.53-yr and/or 6.32-yr pe- 
riodic components. The 5-yr component was also detected 
in CLEANest spectra at 8.0 and 14.5 GHz. Fan et al. [15] 
only detected this 5-yr periodic component from the DCDFT 
spectrum at 4.8-GHz, whereas this component appears at all 
three frequencies in our CLEANest periodograms. The de- 
rived amplitude of the 5-yr periodic component is compara- 
ble to that of the 10-yr component. Thus we regarded it as 
a possible period and subtracted it from the spectrum. Af- 
ter subtracting the three strong periodic components, the re- 
maining peaks in the residual spectra were lower than half 
of the maximum of the original DCDFT spectrum. We fur- 
ther identified two components from the 4.8-GHz spectrum, 
at frequencies of 0.2834 yr -1 , and 0.3518 yr" 1 , correspond- 
ing to periods of 3.52 and 2.84 yr, respectively. The 8.0-GHz 
CLEANest periodogram also presented the 3. 1 1-yr and 2.83- 
yr periodic components. The peaks in the residual spectra 
(bottom panels in Figures [4] and |5]l became less convincing, 
and then the CLEANest stopped. The same procedure was 
applied to the 14.5-GHz data, and we got three periodic com- 
ponents at about 4.20, 3.44 and 2.84 yr. The CLEANed peri- 
odic components are tabulated in Table Q] 



3 Period detection in NRAO 530 



In Section 2 we have demonstrated the powerfulness of the 
CLEANest in fighting against the sidelobes and fake high- 
order harmonics. In this section, we would employ the 
CLEANest algorithm to search for periodicities from the ra- 
dio light curves of NRAO 530. By following the procedure 
described in Section 2, the DCDFT spectra were first derived 
from the observing data. They are displayed in the top panel 
of Figures [4] to [6] The DCDFT spectra show two peaks with 
equal amplitudes at 8.0 and 14.5 GHz (Figs. [5]and[6]), corre- 
sponding to frequencies of ~0.1052 yr"" 1 and ~0.1582 yr -1 . 
The peaks in the 4.8-GHz DCDFT spectrum (Fig. [4]i are 
distributed in a rather broad range from frequency 0.1075 to 
0.1955 yr -1 , indicating a mixture of multiple Fourier com- 
ponents with close frequencies. There are also two high- 
frequency peaks appearing at 0.81 12 yr" 1 and 1.1 140 yr" 1 in 
Fig. H] respectively. The two strongest peaks at frequencies 
of 0.1052 yr" 1 (period = 9.53 yr) and 0.1582 yr" 1 (period 
= 6.32 yr) have already been detected in previous research 
with the same DCDFT method [15], and later confirmed with 
the PDM method [22]. Therefore we used these two periods 
which are confirmed as initial trial periodic components and 
first subtracted them from the DCDFT spectra. In the result- 
ing 4.8-GHz residual spectrum, a Fourier component at fre- 
quency ~0.2000 yr" 1 became the strongest. The two peaks 
shown at 0.8112 yr" 1 and 1.1140 yr" 1 in the DCDFT spec- 
trum became less significant, suggesting they are probably 



In order to verify the detected periodicities, the CLEANed 
Fourier components were used to model the original datasets 
and generate reconstructed light curves of NRAO 530 which 
were compared with the observed light curves in Figure [7] 
The reliability of the detected periodicities is manifested by 
the goodness of the fits of two sets of light curves. At all three 
frequencies, the reconstructed light curves basically recover 
the observed ones from epochs 1980 to 2000. The fits later 
than epoch 2000 and earlier than 1980 show large deviations 
because the data points are sparsely sampled in these time 
spans. Alternatively the deviation between the modeled and 
observed light curves arises because the CLEANed compo- 
nents are fitted with constant amplitudes. However, the ap- 
parent magnitudes of the periodic radio outbursts may alter 
with time considering the complex instabilities of the jet flow 
and/or possibly varying Doppler boosting factors. We also 
tried using different numbers of the periodic components, and 
found that the model light curves constructed using all five or 
six CLEANed frequency components give best match with 
the observational light curves, recovering more fine structures 
in the light curves. 
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Figure 4: The DCDFT (top) and CLEANest (bottom) periodograms de- Figure 6: The DCDFT (top) and CLEANest (bottom) periodograms de- 
rived from 4.8 GHz data. The DCDFT spectrum shows a broad frequency rived for the 14.5 GHz data, 
peak from 0.1075 to 0.1955 yr , indicating the mixture of multiple periods. 
The thick vertical lines in the CLEANest periodogram (bottom) mark the 
subtracted periodic components. 
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Figure 5: The DCDFT (top) and CLEANest (bottom) periodograms de- Figure 7: The reconstructed light curves (thick lines) overlaid on observed 
rived for the 8.0 GHz data. light curves. 
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Table 1: The CLEANest frequency components ofNRAO530 



v obs (GHz) 


# 


Freq. (yr ') 


P(yr) 


Amp. 




1 


0.105263158 


9.50 


1.7536 




2 


0.158227848 


6.32 


1.8280 


14.5 


3 


0.200000000 


5.00 


1.4151 




4 


0.237572585 


4.20 


0.5523 




5 


0.290366492 


3.44 


0.4338 




6 


0.351959385 


2.84 


0.5585 




1 


0.100000000 


10.00 


1.0530 




2 


0.157728707 


6.34 


1.3864 


8.0 


3 


0.195936372 


5.10 


0.9490 




5 


0.321335651 


3.11 


0.6820 




6 


0.352685470 


2.83 


0.6220 




1 


0.107507666 


9.30 


0.7308 




2 


0.146601363 


6.82 


0.7814 


4.8 


3 


0.195468484 


5.11 


0.5677 




5 


0.283429302 


3.52 


0.3781 




6 


0.351843271 


2.84 


0.3405 



In summary, all together five or six possible periods have 
been detected from the NRAO 530 radio light curves by us- 
ing the CLEANest approach. The two strongest periods of 
~9.5-yr and ~6.3-yr have also been reported by other au- 
thors using the DCDFT and PDM techniques [15, 22]. A 
5. 0-yr periodicity is detected in the CLEANest periodograms 
at all three observing frequencies, identified as a possible pe- 
riod. In addition, there are three possible periods which are 
newly found in the present work. They appear at a higher 
frequency regime, corresponding to characteristic periods of 
4.2, 3.4, and 2.8 yr. A possible reason for the new detection 
of three higher-frequency periodic components is that the cur- 
rent analysis uses more data points across a longer time span, 
which results in lower noise in the residual spectra. The rela- 
tively weaker peaks have more chance to be detected. More- 
over these weaker periodic components might not be persis- 
tent across the whole time span. When DCDFT carries out 
a global modeling of the light curves with a linear combina- 
tion of trial functions, any short-lived or intermittent periodic 
components would probably be smoothed in the global power 
spectrum. In the previous search for periodicity in NRAO 
530 [15] the weaker periodic components are probably hid- 
den among the relatively higher noise due to shorter time cov- 
erage and less data points in their study. 

4 Discussion and Conclusion 

A robust Fourier-based algorithm CLEANest was used to 
search for periodicities from the radio light curves of the 



blazar NRAO 530. Five periodic components were detected 
at 4.8 and 8.0 GHz, and six components at 14.5 GHz. The ba- 
sic properties of the detected periodicities can be summarized 
as : 

• multiplicity of the periodicities. The present work con- 
firms that the radio flux densities of NRAO 530 show 
periodic oscillations and there are multiple underly- 
ing periods. Two periods on time scale of 9.5 and 
6.3 yr are most prominent, and they together domi- 
nate the power spectrum. These two periodicities have 
been identified at all frequency bands and with vari- 
ous methods. Monte Carlo tests of the statistical con- 
fidence of these two periods gave a high confidence 
level (Lu et al. 2012). In addition, there are four other 
weaker periodic components, 5.1 yr, 4.2 yr, 3.4 yr and 
2.8 yr. Quantitative statistical confidence estimates of 
these weaker and shorter periods are necessary to clar- 
ify their reliability. A possible reason for the weakness 
of these shorter periods is that they are associated with 
short-lived or intermittent periodic components. Their 
integrated powers over the whole time span are rela- 
tively lower than those of persistent periodic compo- 
nents. Fourier transform carries out a global modeling 
of the entire time series, which tends to wash out the 
(sharp) peaks of short-lived components in the power 
spectrum. These newly detected shorter periods can 
not be the high-order harmonics of the strongest 9.5- 
yr and 6. 3 -yr periods: in each iteration of CLEAN, the 
periodic components identified from the previous iter- 
ation would be subtracted from the DCDFT spectrum, 
while at the same time, the sidelobes and high-order 
harmonic frequency components are removed as well. 
This effect is nicely manifested in the tests of simulated 
time series in Section 2 and Figure|2] 

• mixture of persistent and short-lived periodicities. 
As discussed above, short-lived periodic components 
might be submerged by dominant persistent compo- 
nents. Therefore it would be necessary to seek more 
sophisticated algorithms which have high resolution in 
both time and frequency domains to identify the tran- 
sient periodicities. In fact, in a complementary period- 
icity analysis using the wavelet transform, which is su- 
perior to Fourier-transform methods in the localization 
properties in both time and frequency domains, indeed 
reveals that the 9.5- and 6.3-yr periods are persistent 
across the observing session, while other periodic com- 
ponents only appear in certain time ranges. Mixture 
of multiple periodicities, and the interactions of these 
periodicities increase the difficulty discerning discrete 
individual components. 

• possible harmonic relations of the characteristic fre- 
quencies? The present and previous work provides 
strong evidence for the existence of multiple period- 
icities in NRAO 530 (at least the strongest ones, 9.5 
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yr, 6.3 yr and 5.0 yr are rather significant). If all these 
periodicities are real, an interesting question should be 
answered: do they have intrinsic relationship? Looking 
at the parameters of these periodicities listed in Table 1, 
we find the characteristic frequencies have nearly inte- 
gral ratios. From PI to P6, the ratios are approximately 
2 : 3 : 4 : 5 : 6 : 7. If the relationship of integral ratios 
is confirmed in the future, it would be robust evidence 
that these periodicities have a common origin in phys- 
ical nature. 



The radio emission of blazars is generally thought to arise 
from relativistic jets. A straightforward interpretation of pe- 
riodic radio variability is attributed to the jet precession. The 
periodic outburst is related to periodic ejection of new jet 
components. In this respect, one would expect to see correla- 
tion between the jet ejection and radio outburst. The jet pre- 
cession could be triggered by the orbital motion of a binary 
BH system [26], or a warping accretion disk which results 
from the Lense-Thirring torque due to Bardeen-Petterson ef- 
fect ([27] and references therein). The inferred time scales 
in the above processes vary from hours to years, depend- 
ing on the mass and spin of the black hole as well as the 
orbital properties (e.g., orbital period for binary BH model, 
or transition radius between the inner flat and outer tilted 
disk for warping disk model). However the precession model 
encounters difficulty in the interpretation of the multiplicity 
and possible harmonic relationship of the detected periodic- 
ities. Alternatively, the periodic oscillations of flux densities 
may result from the hydrodynamic instabilities of the jet flow. 
One commonly-seen instability in relativistic jets is Kelvin- 
Helmoholtz (K-H) instabilities [28] that can be triggered by 
either regular perturbations to the jet base (e.g. precession of 
the accretion disk), or random perturbations (e.g., jet-cloud 
interactions). Under certain circumstances, initial perturba- 
tions may grow and propagate along the jet flow as waves 
with various modes such as pinch (n=0), helical (n=l), el- 
liptical (n=2). The hydrodynamic instabilities may shape the 
jet morphology, and also result in periodic flux variations. 
Lobanov & Zensus successfully modeled the emission struc- 
ture of the 3C 273 jet with a combination of helical and ellip- 
tical surface (m=0) and internal (m>l) normal modes [29]. 
Similar approaches were also adopted to study the M87 jet 
[30] and 3C 120 jet [31]. The instability model naturally ex- 
plains the multiplicity and harmonic relationship of the peri- 
odicities in NRAO 530. For instance, the propagating shocks 
in the helically twisted jets periodically pass through the light 
of sight, causing the periodic enhancement of the flux den- 
sity owing to Doppler boosting effect. The observed bright- 
ness distribution of the jet represents a combination of multi- 
ple instability modes. Further monitoring of the jet structure 
variation and temporal variability would be necessary to dis- 
criminate the different models. 
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